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Abstract: Windblown sand transport is a leading factor in the geophysical evolution of arid and semi-arid 
regions. The evolution speed is usually indicated by the sand transport rate that is a function of launch 
velocity of sand particle, which has been investigated by the experimental measurement and numerical 
simulation. However, the obtained results in literatures are inconsistent. Some researchers have discovered 
a relation between average launch velocity and wind shear velocity, while some other researchers have 
suggested that average launch velocity is independent of wind shear velocity. The inconsistence of launch 
velocity leads to a controversy in the scaling law of the sand transport rate in the windblown case. On the 
contrary, in subaqueous case, the scaling law of the sand transport rate has been widely accepted as a cubic 
function of fluid shear velocity. In order to explain the debates surrounding the windblown case and the 
difference between windblown and subaquatic cases, this study reinvestigates the scaling law of the vertical 
launch velocity of windblown transported sand particles by using a dimensional analysis in consideration 
of the compatibility of the characteristic time of sand particle motion and that of air flow. Then a wind 
tunnel experiment is conducted to confirm the revisited scaling law, where the sand particle motion pictures 
are recorded by a high-speed camera and then the launch velocity is solved by the particle tracking 
velocimetry. By incorporating the results of dimensional analysis and wind tunnel experiment, it can be 
concluded that, the ratio of saltons number to reptons number determines the scaling law of sand particle 
launch velocity and that of sand transport rate, and using this ratio is able to explain the discrepancies 
among the classical models of steady sand transport. Moreover, the resulting scaling law can explain the 
sand sieving phenomenon: a greater fraction of large grains is observed as the distance to the wind tunnel 
entrance becomes larger. 
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1 Introduction 


Windblown sand transport is a leading factor in the geophysical evolution of arid and semi-arid 
regions (Bagnold, 1941; Kok et al., 2012; Burr et al., 2015). The launch velocity of sand particle is 
important to quantify the windblown sand transport, to which the experimental measurement is one 
of the main researching approaches. However, the obtained experimental results are inconsistent. 
Some researchers have discovered a relation between average launch velocity and wind shear 
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velocity (Owen, 1964; Bagnold, 1973; Anderson and Hallet, 1986; McEwan and Willetts, 1993), 
while some other researchers have suggested that average launch velocity is independent of wind 
shear velocity (Namikas, 2003; Rasmussen and Sørensen, 2008; Creyssels et al., 2009; Martin and 
Kok, 2017). Numerical simulations conducted by Feng et al. (2009) demonstrated that over an 
erodible sand bed, average launch velocity is independent of wind shear velocity for low wind 
velocities, but increases with wind shear velocity for high wind velocities. Ho et al. (2011) fed sand 
particles at the upwind end of a wind tunnel and experimentally studied the influence of sand bed 
configuration on sand motion. Their results demonstrated that average launch velocity is 
independent of wind shear velocity over erodible sand beds; however, average launch velocity 
depends on wind shear velocity over non-erodible sand beds. The debate on launch velocity has 
been extended to the dependence of steady sand transport rate on wind shear velocity (Martin and 
Kok, 2017), which has been suggested to be cubic (Bagnold, 1941; Kind, 1976), quadratic 
(Creyssels et al., 2009; Martin and Kok, 2017) or even variational (Sauermann et al., 2001; Duran 
et al., 2011; Ho et al., 2011). In contrast to the windblown case, the dependence of steady sand 
transport rate on fluid shear velocity has been determined to be cubic in the subaquatic case (Martin 
and Kok, 2017; Wu et al., 2018). 

The horizontal component of launch velocity is expected to scale with the square root of the 
product of gravitational acceleration and average sand particle diameter according to the balance 
among the forces acting on sand particles. So, it is the vertical component of launch velocity that 
must be responsible for the controversies among the models discussed above. The purpose of this 
paper is to provide a unified interpretation for the controversies based on the scaling law for the 
vertical component of sand launch velocity. The present work is divided to three parts. Firstly, a 
simple dimensional analysis is carried out to revisit the scaling law for vertical launch velocity and 
steady sand transport rate. Secondly, a wind tunnel experiment on windblown sand flow is 
conducted to test the scaling law. Thirdly, an analysis based on the revisited scaling law is proposed 
to explain the controversies among the classical models. 


2 Scaling law revisited by using a dimensional analysis 


The definition of the focus point has to be revisited in the first place. As shown in Figure 1, the focus 
point is the locus where wind velocity is independent of wind shear (friction) velocity (us) (Bagnold, 
1941; Owen, 1964). Sand motion below the focus point is almost independent of wind strength and 
that above the focus point is affected by wind strength. Therefore, it is reasonable to distinguish 
reptation and saltation based on the observed particle trajectory height. Here sand particles that never 
move above the focus point are called reptons and those that can go above the focus point are called 
saltons (Bagnold, 1941; Ho et al., 2014). In subaqueous cases, most sand particles are in suspension, 
which literally move above the focus point. In the present paper, such suspending particles are also 
named saltons for the sake of a unified framework. These definitions are considered the resultant 
expressions of the original definitions, in which the saltons are those driven by wind and impact on 
the sand bed, and reptons are those splashed by saltons (Antdrotti et al., 2002). 
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Fig. 1 Transport layer of windblown sand particles. Ar and hs are the characteristic hop heights (m) of reptons and 
saltons, respectively. 
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Suppose the characteristic times of the sand motion and macroscopic structure of the air flow 
must have the same order of magnitude; otherwise, no interaction would exist between the two 
motions. This is similar to the correlation between the mean flow and the largest eddy in the 
turbulence (Tennekes and Lumley, 1972). As a result, we have: 


v/ g=c,(OU /éz)", (1) 


where v is the vertical launch velocity of sand particles (m/s); g is the acceleration of gravity (m/s); 
Cy is a coefficient of order unity; U is the horizontal wind velocity (m/s); and z is the vertical 
coordinate (m). v/g represents the characteristic time of sand particle motion (s), (OU/dz)! is the 
reciprocal of the velocity gradient and accounts for the characteristic time of the air flow that 
interacts with sand particles. 

The wind velocity at the focus point is proportional to (gd)"”, and it is independent of us (friction 
velocity) below the focus point (Ho et al., 2014). Let's define that the subscripts 's' and 'r' refer to 
the parameters of saltons and reptons, respectively. Then in the reptation layer, the characteristic 
wind velocity U, scales with (gd)'”, and the characteristic length is the hop height of the reptons, 
h, which scales with v} /(2g). Following Equation 1, the vertical launch velocity (vr) of reptons is: 


v, ~ g(h,/U,) ~ v1 (2(gd)""”) . (2) 


wom 
~ 


So one can obtain v:=2c:(gd)!?, where "~" means that the scaling factor between the two sides is 


within a magnitude of five times (Tennekes and Lumley, 1972). That is to say, cris a coefficient 
that ranges from 0.2 to 5.0. 

In the saltation layer, as shown in Figure 1, the characteristic length is the hop height of the 
saltons As, which scales with v,/(2g). The characteristic wind velocity U, scales with the effective 
wind shear velocity (umur), where ug is the threshold wind shear velocity above which sand 
particles can be set in motion. According to Equation 1, the vs of the saltons is: 


v, ~ g(h, 1U) ~ v? 1 (2(u; —Uug)) - (3) 


So vs=2cs(ur-un). Similarly, cs is a coefficient ranging from 0.2 to 5.0. As a result, the hop heights 
for reptons and saltons can be expressed as: 


h, =2c?d 
h, =2c? (u; —up) 1g 
The hop length is /Æum(2v/g), where um is the average horizontal velocity (m/s) of sand 
particles. Although the characteristic horizontal velocity of sand is unknown, moving sand is 
accelerated by wind in the horizontal direction. Therefore, um is roughly proportional to the 
characteristic wind velocity. Then one obtains: 
l cU (2v,/g)œ£d 
l, c U (2v, / g) x (uy =u}? 18 


where /, and /; are the hop lengths (m) of reptons and saltons, respectively; and œ means that the 
left side is proportional to the right side, and the range of the scaling factor is uncertain. Comparison 
of Equations 4 and 5 shows that the ratio of hop length to hop height is constant for reptons and 
saltons, which is consistent with experimental results (Zhang et al., 2007; Ho et al., 2011, 2014). 
On the basis of momentum conservation, the steady sand transport rate can be obtained as 
(Sørensen, 2004): 


(4) 


(5) 


q=l1,(t-1,)/ (u, =u), (6) 


where Zm is the average hop length (m); t is the shear stress (t=pu} with p=1.25 kg/m? being the 
air density); t is the threshold shear stress (t1=puj.); ut and u, are the average horizontal velocities 
of sand particles taking off from and landing on the sand bed, respectively; and (u-u+) is roughly 
proportional to the characteristic wind velocity as (u\-u;)*<(gd)' if the majority of the moving 
particles are reptons, and as (u-u1)°= (usut) if the majority of the moving particles are saltons. 
Combination with Equation 5 yields: 
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q, < pu; —u;,)(d/g)"”, (7) 
q, © PU; —Uy (Up —Ug)/ g- (8) 


3 Wind tunnel experiment 


Equations 4-8 are derived by supposing that vertical launch velocity scales with (gd)! for reptons 
and scales with (usur) for saltons. A wind tunnel experiment is introduced to testify this 
conclusion. 

The experiment is conducted at the Laboratory of Windblown Sand Dynamics, Xi'an Jiaotong 
University, China. As shown in Figure 2a, such laser based configuration has been widely used in 
the measurement of windblown sand transport (Zhang et al., 2007, 2008; Wang et al., 2009; Yang 
et al., 2009; Kang, 2012; Zhang et al., 2014; O'Brien and McKenna Neuman, 2016, 2018). The test 
section of the wind tunnel has a size of 6.75 m*x0.72 m’x0.60 m’, with x denoting the mainstream 
direction, y the transverse direction and z the vertical direction. The wind is sent out by two abreast 
air blowers. In order to fully mix the two air jets and to smooth the turbulence caused by air blowers, 
a distance of 5.0 m, which is paved with roughness elements, is set between the wind tunnel 
entrance and the upstream edge of sand bed. Thus the sand bed with a size of 1.0 m*x0.2 m’x0.01 
m’ is placed 5.0 m downstream of the wind tunnel entrance. A high-speed digital camera 
(FASTCAM-APX 120K, Photron, Japan) is placed at the side of wind tunnel, with its view line 
perpendicular to the side wall plane. The frame rate of the digital camera is 2000 Hz, the exposure 
time is 1/3000 s, and the resolution is 1024 pixelx572 pixel. The observation region is 138.4 
mm*x69.2 mm’, which is illuminated by a sheet laser (Mini-YAG, 100 mJ). The diameter range of 
the three sand samples in the present experiment are 100-125 um with a density of 2702.8 kg/m’, 
200-300 um with a density of 2650.0 kg/m? and 300-500 um with a density of 2650.0 kg/m?. One 
of the recorded particle images in the present experiment is shown in Figure 2b. The moderate sand 
particle concentration (especially above the focus point) allows for the Lagrangian particle 
detection to catch the particle trajectory. Firstly, the dynamic threshold scheme (Ohmi and Li, 2000; 
Zhang et al., 2016) is used to realize the image segmentation and the recognition of moving 
particles. On a grey-scale particle image, a small window is built and shifted, and a local grey-scale 
threshold in this window is continuously computed and applied to distinguish the particle image of 
multiple grey-scale levels from background (Fig. 2c). This operation is also called the 
reconstruction of particle image. Secondly, after the operation of Voronoi diagram on a 
reconstructed particle image sequence, each particle across the sequence corresponds to a unique 
Voronoi polygon, so each particle can be tracked across the frame based on the geometric stability 
of such polygon (Zhang et al., 2015, 2016). As a helper method, a central-differential schemed 
static background elimination method called the "non-zero subtraction algorithm" is applied to 
extract the moving sand particles near the sand bed (Zhang et al., 2016). Nevertheless, those 
touching closely the sand bed are difficult to extract due to the over exposed reflection at the sand 
bed. The result is satisfactory as shown in Figure 2d. 

The streamwise mean velocity profile U=U(z) is measured by an array of hot-wire probes 
(55R49, Dantec Inc., Denmark), which are the nikel-plating type and located downstream of the 
observation region. Before the measurement, the sand bed was fixed by spraying water. The 
reference wind velocity Uo is measured 0.26 m above the sand bed surface. As shown in Figure 3, 
the profile is well-fitted by the logarithmic law. 

The wind shear velocities are calculated from the wind velocity profiles, which are 0.23 and 0.30 
m/s for the two cases of 100-125 um sand, 0.40 and 0.46 m/s for the two cases of 200-300 um 
sand, and 0.41 and 0.51 m/s for the two cases of 300-500 um sand, respectively. The particle 
Reynolds number Re are 1.7136 and 2.2351 for 100—125 um sand, 6.6225 and 7.6159 for 200-300 
um sand, and 10.8609 and 13.5099 for 300-500 um sand. The characteristic length and velocity 
are sand diameter and wind shear velocity. In another independent set of tests, the threshold velocity 
of sand motion for the above three sand samples is determined by using the technique shown in 
Figure 2d, and the corresponding threshold wind velocity profile is measured and used to calculate 
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ur. The results are 0.15, 0.23 and 0.29 m/s for the samples corresponding to the diameter ranges of 
100-125, 200-300 and 300-500 um, respectively, agreeing well with the classical model 
up=0.1[gd(ps-p)/p]'”, where p=1.25 kg/m? is the air density. 


Fig. 2 Measurement of windblown sand flow in wind tunnel. (a), physical photo of the experimental setup; (b), 
one of the recorded particle image and its partially enlarged display; (c), recognition of the connected domains by 
dynamic threshold scheme. The dark blue represents the background, and the different colors represent the different 
levels of image segmentation; and (d), trajectory of sand particles. 
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Fig. 3 Streamwise wind velocity (U) profiles. The dashed lines refer to the linear fittings of the relationship of U 
and In(z). 


The maximum height of the sand particles in each image is approximately the hop height of the 
saltons. Therefore, the average value of the maximum heights for all particle trajectories is defined 
as hs. The average height for all moving sand particles is defined as Am. These two kinds of height 
of moving sand particles versus the effective wind shear velocity are shown in Figure 4. To 
determine if the windblown sand flow has reached the steady state, each image is divided into two 
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equal parts along the flow direction, and the change rate (R) of moving sand particles is calculated 
as R=(Nu-Na)/(Nut+Na), where Nu and Na are the numbers of moving sand particles at the upstream 
and downstream parts of the region covered by the recorded image, respectively. When most of R 
is less than 0, it means that more and more sand particles are blown up as wind moves toward 
downstream, i.e., the sand flow does not reach the steady state. However, when R fluctuates around 
0, it means that sand particles carried by wind reaches the dynamic equilibrium, i.e., the sand flow 
is in the steady state. In this case, the fluctuation of R results from the intermittent nature of sand 
flow (Platzer and Francois, 2011). The result of R is shown in Figure 5. As one can see, for each 
sand sample, the sand flow reaches the steady state only at the relatively low wind shear velocity. 
According to Figure 4, however, the hop height of saltons and the average height of all moving 
sand particles are related to the wind shear velocity only, thusly disregarding whether the sand flow 
is steady. Therefore, the unsaturated state does not affect the scaling law of launch velocity of sand 
particles. 

120 
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Fig. 4 Hop height of saltons and average height of moving sand particles versus the effective wind shear velocity. 
The height is normalized by mean diameter d, and the effective wind shear velocity squared is normalized by gd. 
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Fig. 5 Change rate (R) of the number of moving sand particles along the flow direction versus time. (a) and (b) 
are for the 100-125 um sand with u=0.23 m/s and u=0.30 m/s, respectively. (c) and (d) are for the 200-300 um 
sand with u=0.40 m/s and u=0.46 m/s, respectively. (e) and (f) are for the 300—500 um sand with u=0.41 m/s and 
ur=0.51 m/s, respectively. 
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According to the results shown in Figure 4, the relationship between the dimensionless hop 
height and effective shear velocity can be fitted linearly as 


h./d=4.37(u, —u,)° /(gd)—0.20 
h,, 1d =0.34(u, -u,)? I (gd) +3.74 

By comparing Equations 4 and 9, we can obtain c;=1.48. The average height is 
h,,/d=Ph,,/d+(Q-P)h,,/d, (10) 


where the average height of saltons Ams equals to 0.5hs because the heights of saltons in each image 
are assumed to be distributed evenly in the range of 0—As. Similarly, the average height of reptons 
hmr is 0.5%; P is the percentage of saltons among all moving sand particles and remains 
approximately constant because the explored range of uput, namely, 1.41—2.00, is narrow. Am/d and 
h,/d are linear functions of (up-wp)?/(gd). According to Equation 10, the experimental result of /,/d 
can be also fitted with a linear function, which is assumed as: 


h,/ d = Au, —u,) /(gd)+B, (11) 


where A and B are unknown coefficients. By substituting Equations 9 and 11 into Equation 10, one 
can obtain the following relationship: 


(0.68 —4.37P — A(1 — P))(u, u)? /(gd)+7.48+0.20P -B(1-P)=0. (12) 
This equation should hold for all values of (wp-wn)*/(gd), that is: 
0.68 —4.37P— A(1- P) =0 
7.48+0.20P-BU-P)=0 


The sand motion below the focus point is independent of wind shear velocity (Ho et al., 2014). 
Hence, A=0 according to Equation 11. Then, P=0.16 and B=8.94. The fitted line for the 
experimental result of A/d is h,/d=8.94, which yields c,=2.11 according to Equation 4. Therefore, 
v=4.22(gd)!? and v;s=2.96(ur-un). As shown in Figure 4 and Equation 9, the slope of hm/d versus 
the effective shear velocity squared is nearly zero and about ten times smaller than that of h;/d; 
meanwhile, the slope of /,/d is theoretically zero. That is to say, reptons contribute to the most of 
the average hop height. This is consistent with the fact that sand particle concentration decreases 
drastically (such as exponentially) as the height increases (Dong et al., 2003; Kang, 2012). By 
using a sand sample with a mean diameter of 230 um, Ho et al. (2011) found that the average 
vertical launch velocity is 10.64(gd)'? over an erodible sand bed, whereas the average vertical 
launch velocity is 5.8(u}-ur) over a non-erodible sand bed. As discussed in Section 4, the majority 
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Fig. 6 Trajectories of sand particles on the flat sand bed. Compared to the grain in the present experiment, the 
grain in Ho et al. (2011) will experience an extra acceleration process before the first launch process. 
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of the moving sand particles are reptons for the erodible sand bed and are saltons for the non- 
erodible sand bed. Therefore, their scaling law is qualitatively the same as that of the present study, 
while their coefficients cr (5.32) and cs (2.90) are larger than the corresponding values in the present 
experiment. This is because the sand particles in their experiment was supplied at the upwind end 
of the wind tunnel, which resulted in an extra acceleration process before the first launch process 
(Fig. 6). As a result, the launch velocities as well as the hop lengths are larger than those in the 
present experiment. Above all, the scaling law in Equation 4 are consistent with the results extracted 
from Figure 4, the revisited scaling law are thusly considered reasonable. 


4 Explanation of the discrepancies among the classical models 


According to Equations 7 and 8, the dependence of q on up varies from quadratic to cubic as the 
ratio for particle number of saltons to reptons increases. These scaling law is used to explain the 
discrepancies among the classical models of q shown in Table 1. In a word, the scaling law of q 
depends on the ratio of saltons number to reptons number. For a small uf over erodible sand beds, 
reptons account for the majority of moving particles, leading to a quadratic dependence of q on uf 
according to Equation 7, which agrees with the results of Creyssels et al. (2009) and Ho et al. (2011) 
for erodible sand beds, Duran et al. (2011) for small uf, and Martin and Kok (2017) for a large ur 
over erodible sand beds or for non-erodible sand beds where the height of the focus point is only 
d/30, saltons are the majority. Therefore, the dependence of q on uf should be cubic according to 
Equation 8, thereby supporting the results of Duran et al. (2011) for large uf and Ho et al. (2011) 
for non-erodible sand beds. The dependence of q on us, which varies from quadratic to cubic 
(Sauermann et al. 2001), is caused by the increase in the ratio of saltons to reptons. The cubic 
scaling law for q in Bagnold (1941) and Kind (1976) was obtained by assuming that the vertical 
launch velocity scales uniformly with uf. However, this assumption is unconvincing because the 
vertical launch velocity of reptons is independent of ur. 

The cubic scaling law for g in subaquatic cases can be understood in the same way. In windblown 
sediment transport, the wind velocity profile is logarithmic above the focus point and is independent 
of urs below the focus point (Ho et al., 2014). By contrast, the fluid velocity profile is logarithmic 
throughout the boundary layer (Lajeunesse et al., 2010). As a result, water-driven sand flows have 
no focus point, that is, the height of the focus point is zero. Therefore, all moving sand particles 
can be viewed as saltons, leading to a cubic scaling law for q. 


Table 1 Models for the steady sand transport rate 


Reference Expression ug Ug 
Bagnold (1941) qœ (dldy)'?pui z 
Kind (1976) q= puku: -ur )/ g - 
Creyssels et al. (2009) q=<p(u;—u;)(d/g)\” <3.4 
Martin and Kok (2017) qœ plui -ur )(d/g)"? <2.2 
Sauermann et al. (2001) q%ur, nE(2, 3) - 


œ p(u2—u;,)(d/g)'” for erodible bed z 
Ho et al. (2011) q5 plur )die) l 
q&<p(ur—Uug)( ut—un)/g for non-erodible bed - 


œ puç (ur —ui,)/g at large ur - 
Duran et al. (2011) T p i i i ' 
q&p(u;—ug’)(d/g)"~ at small ur - 
Note: do=250 um is the reference sand diameter. Sauermann et al. (2001) suggested that q is proportional to ur, where n increases from 2 
to 3 as the wind shear velocity increases. -, no data available. 


5 Conclusions 


By incorporating a dimensional analysis and a wind tunnel experiment, the present work shows 
that the ratio of saltons number to reptons number determines the scaling law of sand particle launch 
velocity and steady sand transport rate. When saltons account for the majority of moving sand 
particles, the average launch velocity depends on wind shear velocity, and the steady sand transport 
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rate scales with the cubic of wind shear velocity; when reptons take the majority, the average launch 
velocity is independent of wind shear velocity, and the steady sand transport rate scales with the 
quadratic of wind shear velocity. This conclusion provides a succinct framework in which the 
differences among the previous experimental results can be interpreted. Furthermore, for a 
moderate shear velocity, reptons take the majority of moving sand grains for aeolian transport. Then 
according to the formula of hop length /,-<d, one can infer that in aeolian transport the larger the 
sand size is, the farther the sand jumps. This might be the reason for the sand sieving phenomenon: 
a greater fraction of large grains is observed as the distance to the wind tunnel entrance increases 
(Duran et al., 2011). The next step is to investigate a more effective way of monitoring the temporal- 
variant ratio of saltons to reptons under a wide range of wind velocity, so that the windblown sand 
transport can be correctly quantified under various circumstances, and the critical velocity at which 
the scaling law switches can be distinguished. 
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